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Abstract 

We study the propagation of leptons in hot and dense media and 
estimate the background corrections to the purely leptonic processes 
in the early universe and in the stellar cores. We show that the first 
order radiative corrections to purely leptonic processes in the real-time 
formalism can easily be evaluated without going into the resummation 
of thermal loops in these media. It is also possible to express the hot 
and dense corrections to the leptonic scattering cross-sections in terms 
of the corresponding corrections to the decay rates of the (heavy) 
intermediate particles in the real-time formalism upto the one loop 
level. 

1 Introduction 

It is now well known that the stellar cores can have very high temperatures 
and densities [1]. It therefore looks more reasonable to incorporate the finite 
temperature and density (FTD) effects in the study of stellar media. The 
physics of superdense collapsing stars is expected to be better understood 
incorporating the statistical corrections due to the hot and dense background 
in the stellar cores. The superdense collapsing stars like supernovae have core 
temperatures around 30-70 Mev and the chemical potential of the electrons 
\i e ~ 250 — 300 Mev [2], which is sufficiently below the masses of W and Z. So 
the important contributions from the background can be obtained from the 
hot and dense leptons in the background heat-bath at thermal equilibrium. 
The standard big bang model of cosmology [3] also predicts the creation and 
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annihilation of leptons in the very early universe. These processes were taking 
place during the first few seconds of the creation of the universe and leptons 
were the first matter particles. The cross sections of these reactions are also 
related to the expansion rate of the universe [3]. The temperatures of the 
universe at that time were high enough to show non-negligible contribution 
of the hot and dense lepton background. The contributions due to the hot 
electroweak massive bosons is not included because the background particle 
energies are not sufficiently high to have real W's and Z's in the background. 
Moreover, the lepton background effects will always dominate over these 
heavy mass particles. 

We discuss here the pure leptonic decays and cross-sections in the Stan- 
dard Electroweak Model (SEWM). We will not discuss background radiative 
corrections to pure QED processes (Bhabha and Moller Scatterings) here, 
because they deserve a detailed independent study of the topic, which is not 
in the scope of this paper.. We mainly evaluate the first order radiative 
corrections to the neutral current as well as the charge current processes in 
hot and dense lepton background. Moreover, the first order radiative correc- 
tions to the inelastic and to the elastic scatterings of leptons in the standard 
model are discussed through the neutral current or the charged current pro- 
cesses taking into account the thermal equilibrium of leptons and photons 
in the background . Strictly speaking, we restrict ourselves to the SEWM 
so the additional corrections to these processes from the extended standard 
model [4,5] are not incorporated in the same statistical background. For 
this purpose we will have to consider the massive neutrino physics in such a 
medium. We restrict ourselves to the temperatures where the light particles 
can only be found in the background heat-bath and the first order QED type 
radiative corrections are significant. We use covariant real time formalism 
[6,7] in which the renormalization of QED at finite temperature and density 
(FTD) up to one loop level [7,8] have been discussed in detail, for all pos- 
sible ranges of T and [i. It has also been checked that the renormalization 
scheme of QED in vacuum works perfectly at FTD and there is no need 
of resummation [9] at the one loop level in QED for real particle systems. 
It is mainly because of the reason that at the sufficiently low temperatures 
the hard thermal loops cannot appear as the background temperature is not 
high enough to generate real hadrons to incorporate the statistical corrections 
from hadronic background. Moreover, the box diagrams with the overlapping 
phases do not contribute [10]. Therefore, the 1-1 component of the propaga- 
tor is important considering the physical particles in the background [7,8]. 
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The renormalization constants of QED have been evaluated upto the two 
loop levels in literature [11]. We don't need to incorporate the overlapping 
phases, because the overlapping phase contribution is automatically approx- 
imated out. These additional thermal contributions can then be estimated 
easily by applying the renormalization techniques of QED. This procedure is 
not similar to Cutkowski Rules [12]. We are not using any cutting of diagram 
as a rule, it is the approximation under which we keep the leading order 
contributions and ignore the rest. We could do it only because the real time 
formalism give thermal/statistical correction term as the separate additive 
terms in all kind of Feynman diagrams. We use the renormalization scheme 
of QED and since the angle of scattering is independent of statistical back- 
ground parameters, we can simply express the cross-sections in the form of 
the decay rates. Using the values of these renormalization constants evalu- 
ated in the real time formalism we can easily evaluate the first order radiative 
corrections to the decay rates and the scattering cross-sections. The first or- 
der radiative corrections to the decay rates of H° — > fl [7,8] and Z° — > ii 
[13] have already been calculated for temperatures above the chemical poten- 
tials. We generalize these techniques to evaluate the background corrections 
to the leptonic decay processes of real particles and give a general discussion 
of the QED type corrections to the electroweak vertices at FTD. Section 2 
comprises the calculational scheme, whereas, we give the analytical forms 
of the radiative corrections in different limits of temperatures and chemical 
potentials in section 3. 

We have made use of all the existing results of the QED renormalization 
constants and filled up the gaps wherever they were. We had made use of 
the renormalization constants to estimate the decay rates of the relevant 
particles and then related the lepton scatterings to these decay rates. All the 
calculations are done in the standard electroweak model. Some of the possible 
applications of these results in astrophysics and cosmology are discussed in 
Section 4. 

2 Calculational Scheme 

The statistical background corrections to the physical processes are 
calculated using the Feynman rules of vacuum theories. However, the particle 
propagators are modified in hot and dense medium using the distributions of 
the background particles. The boson propagator in the real-time formalism 
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can be written as 
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k 2 + ie 

and the fermion propagator as, 
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with the Bose-Einstein distribution, 
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for massless vector bosons. This automatically leads to the vanishing mass 
density and hence the chemical potential [14]. The Fermi-Dirac distribution 
is given by [15] 
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such that the chemical potential of fermions is always equal and opposite to 
the corresponding antifermions. Since all the particle species are in thermal 
equilibrium, it is convenient to expand the distribution functions in powers 
of (3po, where (3 is the inverse temperature and po is the energy of the cor- 
responding particles. T, n and m are all expressed in the same units. We 
obtain[16] 



np (-l) n e n/ ^ K >) 



(5a) 



n=l 



and 



np (-l)"e- n ^-P°) + 0{n - po) 



(5b) 



n=l 



for particles, and 



n F ^J2 (-l) n e- n ^ +W) ) 



n=l 



(5c) 



4 



and 

oo 

np (-l) n e - nl3 ^ +Po) (5d) 

n=0 

for antiparticles in the corresponding temperature and density regime. We 
do all our calculations by the corresponding substitution of the distribution 
functions from Eqs.(5). We have explicitly checked the validity of the renor- 
malization scheme of QED by reproducing the physically measurable results 
of the first order background corrections to the decay rates by adding all the 
one loop diagrams with those by replacing the vacuum parameters with the 
corresponding renormalized parameters at FTD upto the one loop level. 

3 Calculational Details 

We discuss the renormalization constants of QED in hot and dense medium. 
We mention here the results of the first order radiative corrections to the 
self mass, wavefunction and the charge of electron (QED coupling). We 
have used the renormalization scheme of QED in the calculation of radiative 
corrections to the gauge vertices of electroweak theory. We have explicitely 
demonstrated the order by order cancellation of singularities upto the two 
loop level in QED. We have calculated the QED type of hot and dense back- 
ground corrections to the Z° — > il electroweak vertex to explicitely check 
the renormalizability of the corresponding gauge theories and to evaluate 
the radiative corrections. The renormalization constants of QED have al- 
ready been calculated upto the two-loop level in our scheme[ll]. We will 
be using this renormalization scheme in our calculations to estimate the sta- 
tistical background corrections. We briefly mention the previously calculated 
renormalization constants of QED and tried to put everything together in 
the convenient form and have tried to fill the gaps if there were any. We 
refer to the existing results for the details, wherever it is relevant. Then 
using all these results we go to the calculations of the decay rates. Then 
we discuss, how we can use the radiative corrections to the decay rates to 
evaluate the first order radiative corrections to the crossections. We have to 
rewrite the already existing results sometimes, to develop a uniform way of 
writing everything for the comparison purposes. 
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3.1 Renormalization of QED 



We calculate the QED type radiative corrections here so we use the renor- 
malization of QED and make use of the renormalization constants given in 
literature [4-18] in the relevant ranges of T and fi. We only rewrite them in 
the convenient form. The FTD radiative corrections to the electron mass — 

m 

correspond to the renormalization constant of electron at T > fi and can be 
expressed as [7] 
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and, in the limit T < fi , — appears [17] to be 
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for electron. It is also worth-mentioning that for T > fi the electrons and 
positrons get equal contribution from the background. However when T < fi 
the situation is different. In the classical limit T <C m < fi and A,B,C,and 
D functions vanish for electrons and the chemical potential correspond to 
the Fermi energy. In this case the there will be no antiparticles in the back- 
ground and we don't have to bother so much about the self-mass corrections 
to positrons. However, just for completenes we evaluate the background 
contribution to the positron mass at T < /i as 
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The wavefunction renormalization constants can be obtained from the 
selfmass of the corresponding particle using Ward identities [18]. We have 
obtained the renormalization constant of the e~ and e + wavefunction at T 
> /it, as 
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whereas, in the limit T < fi one obtains different forms of the wavefunction 
renormalization for particle and antiprticle. The e~and e + wavefunction in 
the high chemical potential limit come out to be [17] 
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The positron wavefunction renormalization can easily be evaluated by replac- 
ing the distribution function of electron by that of positron in highly dense 
systems. It comes out to be 
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The electron charge renormalization in QED is determined through the 
self-energy of photons. So, the fermion charge renormalization constant in 
the limit T > /i,it comes out to be [18] 
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The fact that QED is renormalizable in this regime and also there is no 
doubt about the renormalizability of electroweak theory because no extra 
hot infrared singularities enter into the calculations due to the massive in- 
termediate vector bosons. All the background corrections to the electroweak 
processes are suppressed by the square of the electroweak gauge boson masses 
as compared to the lepton masses so we don't have to incorporate them here. 

It has been explicitely checked that the QED renormalization scheme 
works perfectly in the quantum statistical medium. The calculations of the 
leptonic decays, using the renormalization constants at FTD reproduces all 
the existing results obtained by adding all the first order diagrams in QED. 
In the next subsection we discuss the leptonic decay rates in hot and dense 
background in a little bit more detail to see that how can we use these decay 
rates to evaluate the leptonic crossections. 
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3.2 Leptonic Decays 

The radiative corrections to the decay rates of charged leptons can be es- 
timated from the renormalization constants of QED . The decay rate of Z° 
(i.e.e;the Z width) is an important gauge vertex of the electroweak theory. 
We have evaluated it in the limit T> u in the rest frame of Z°. It comes out 
to be equal to [14] 
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and the scalar Higgs decay rate is calculated in the rest frame of Higgs which 
comes out to be equal to [17] 
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mii is the mass of the neutral Higgs of SEWM. We generalize the existing 
results of Ref.[17] to evaluate the Higgs decay rate at high densities at non- 
negligible temperatures. We use the approximation m 2 H » m 2 to ignore 

the higher order terms in . The calculational techniques are same as are 

used in Ref. [17]. This can be considered as a straight forward generalization 
because all the techniques are same. The neutral Higgs decay rate for li>zT 
, then appears to be equal to 
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We use the rest frame of Z° and also the energy of Z° is high enough to 
ignore higher order contribution from the distribution function. We include 
all the one loop Z diagrams including the bremstrahlung. The calculational 
scheme is very similar to the one adopted previously in Ref. [13,17]. We con- 
clude that the leading order contributions due to the highly dense background 
comes out to be some complicated function of several parameters given be- 
low. They are even more complicated than the results of Eqs.(19 & 20). We 
postpone the detail discussion of these results for now. 



The calculational scheme which we have been using throughout can ac- 
tually help us to evaluate the statistical corrections to any leptonic decay 
processes. We, therefore, evaluate the QED type radiative corrections to the 
meson decay in the statistical background. We make use of the renormaliza- 
tion constants of QED upto the one loop level and replace all the relevant 
bare parameters with the corresponding FTD corrected values. 

The thermal corrections to the decay rates of neutral mesons like pions 
and kaons can also be calculated in a similar way. They simply comes out 
to be very similar to the thermal corrections to the Higgs decay rates. This 
happens only when we cannot consider any corrections to the vertices other 
than QED type radiative corrections. It appears to be an over-simplification 
of the mathematical calculations but the results are acceptable in the sta- 
tistical theory where we need to assume real highly relativistic hot particles 
in the background. This is the beauty of thermodynamics that we have to 
be in the systems where we can deal with the real particles on mass shell in 
the heat bath. Now the heavy particles exchanging contributions can easily 
be ignored as compared to the lighter ones based on the expected relativistic 
behavior 
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Fi, F2, and F 3 are very similar to the corresponding coefficients in Eq.17. 
Main difference occurs in the bare decay rate, but the radiative corrections 
(especially the leading order QED type contributions) are fairly similar and 
so are easy to evaluate for leptonic decays of mesons also. We have always 
preferred to use the rest frame of the decaying particles to evaluate these 
contributions. 

However, 4E 2 = vr? meson , whereas, m meson correspond to the mass of 
the decaying meson and £ correspond to the lighter leptons(ni£ -< mmeson). 
For heavy mesons the above relations can also hold for taons. In the lep- 
tonic decays of charged leptons, the vertex corrections terms do not get any 
background contribution. So, we can simply replace rest masses and the 
vacuum couplings by the corresponding physical masses and the background 
renormalized couplings. The radiative decays of the charged mesons should 
however appear like 
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In these kind of three body decay processes, the phase space corrections 
due to the absorption and emission of photon at the electron leg of this 
vertex contribute. We can evaluate Eq.23 for all the interesting ranges of 
temperature and chemical potential using the corresponding expressions of 
n F (p T- /-0 • All the physical parameters, like electron charge mass and the 
wavefunction are needed to be replaced by the corresponding renormalized 
(background corrected) values. These relations are mentioned in the next sec- 
tion in Eqs.24. We can evaluate the radiative corrections to different leptonic 
processes upto the one loop level in almost all the ranges of temperature and 
chemical potential in our frame of calculations. However, leptonic scatterings 
are very important in the early universe when the temperatures were suffi- 
ciently smaller than the electroweak gauge bosons mass and still high enough 
to use the limit T^> m e , for m e to be the electron mass in the electron dis- 
tribution function. We can obviously incorporate the FTD corrections using 
our frame of calculation. Since the high density corrections, where densities 
are much larger than the masses, we have to be careful to ignore FTD cor- 
rections which are not QED type. Also, in those dense systems, we might 
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need to incorporate the magnetic field effects also. So, we postpone the dis- 
cussion of lepton scatterings in the highly dense systems at the moment. In 
the next section, we evaluate the lepton scatterings at high temperature and 
low densities and try to investigate that how the background corrections can 
effect the lepton scattering cross-sections in the early universe and also their 
implications to the big bang cosmology. 

3.3 Scatterings of Leptons in Electroweak Theory 

In the calculation of the cross-sections, one can easily make use of the renor- 
malization constants of QED evaluated in hot and dense media. The leptonic 
scatterings can either take place in QED or in the electroweak theory. QED 
is important for charged particle processes of ££ — > ££ whereas £ correspond 
to the charged leptons only. In the standard model additional leptonic pro- 
cesses of the type — >• £^1/ gor ££ — > z/77.would also be expected. The 
individual lepton number conservation would be required in the standard 
model. A very quick estimate of the background corrections can be made by 
simply replacing the bare lepton masses, charge of electron(QED couplings) 
and the wavefunctions by the background renormalized mass, charge, and 
wavefunctions correspondingly. So, we can directly replace the bare mass of 
lepton be the corresponding renormalized mass of leptons given as 

m e -> m e [l + — ] (24a) 

me 

The wavefunction of charged fermions by the corresponding renormalized 
wavef unction given as 

■$>e -»• Z^e (246) 
and the coupling constant or the charge of electron can be rewritten as 

e -> Z^e (24c) 

The values of all the coupling constants can be used for the given limits 
of temperature and the chemical potential from previous sections. The phase 
space corrections can also be obtained from the cross-sections of the corre- 
sponding processes directly. Just to keep ourself limited to the electroweak 
scatterings of particles in the relevant temperature and density regime, we 
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discuss in a little more detail, the electroweak scattering of leptons in the 
physical conditions of the early universe. 

The purely leptonic scatterings can either take place in QED or in the 
electroweak theory. In QED, the elastic or inelastic scatterings of charged 
leptons, namely ££ — > ££ can take place for £ = £ or £ ^ £. We have ex- 
plicitely checked that the QED renormalization scheme works perfectly well 
in hot and dense medium upto the one loop level. Due to this reason, it 
makes no difference if we calculate each and every first order diagram sepa- 
rately and add them together or calculate the matrix elements directly from 
the renormalization constants we reach the same result. It has also been 

2 

noticed that the radiative corrections are obtained in a power series of -^f 

m x 

where nix correspond to the mass of exchanging heavy particle. We use this 
approximation because the heavy particle exchange in the standard basically 
allows the exchange of , Z° , and the standard model Higgs H°. These 
particles have such a heavy mass values that we can easily ignore the higher 
order terms in the power series of ^f. Due to the presence of this small 
parameter, we have to restrict to the first order radiative corrections. Higher 
order in the coupling constant makes the radiative corrections fairly negligi- 
ble. Moreover, the temperature of our system is small enough to make the 
terms proportional to even more negligible. Starting from the matrix 

m x 

element for ££ — > ££ ,we have to replace the ordinary vacuum fermion and bo- 
son propagators by the corresponding FTD propagators from Eqs.l&2. Now, 
using the properties of Dirac S, we can only get the nonzero contributions 
which are of the same form which we obtain in the decay rates because the 
radiative corrections from one vertex can be considered at a time to be in 
the lowest loop level. .In this situation, the statistical contributions to the 
cross-section can directly be related to the decay rate upto the first order in 
a. We can write it as 

daP(££^£T) n TP(Z° ^ £l) t n TP(H° ££) 



do~r 



+ 2 : 



-> ££ ) T (Z° -> ££) T (H° -> ££) 

At the very low energies, where exchanging particles have very low mo- 
menta, the contribution of all of these diagrams is added up together to 
give 



daP(££ -> ££ ) 4aT 
da (££ -> ££) ~ nm z v 



2m £ 
1 + 

m z 



m 2 z v 2 



£n 



1 + v 2ni£ 
1 — v nizv 



3 + 



m 2 z v 3 



a(m e /3,fj,) + 
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4a 

7T 

4a 

7T 



1 + 



2m\ 
m\v 2 



5mf 
4m\v 2 



+ 



mi 



4m 2 z v 4 



1 + v 



2T 2 



m%v 



ml 



m%v 3 



, l + v 1 
1 — v v 



2m 2 



m\v 2 



b(m e (3, n) + 
c(m e (3, fi) 



+ 



4aT 



nm H v 
2a 



-m e + 



7TV Z 



3v< 



8aT 2 



irv^m 



H 



in 



vln 
1 - 



l + v 

v 



m t - m H - 



V + 22 m; 



1 + V 



in 



1-v 



m 



H 



1 + V 

b(m e (3, a) - 



a(mef3, fj) + 



l + v 



+ 3v 



c(me/3, //) 



(26) 



which gives the leading order contribution at very high temperatures as 



da f} {Zi-*M)„2ot'nT 2 

da (ee -»• if) v 2 



1 



3m| 



+ 



1 



rn 



H 



(27) 



for the photon energy E 2 » mj, and T > /i. 

The radiative corrections to the scattering cross-sections of charge current 
processes can be evaluated in terms of the decay rates of W mainly. 



r (^± -> £±^) 



This contribution is always proportional to 



T 2 



(28) 



-.So, larger contribu- 



tions are obtained for larger temperatures. 

For the limits of T < //, in the relativistic limit, it turns out to be a com- 
paratively complicated expression. We don't want to involve in the detailed 
discussion of these contributions because in the early universe and in the 
stable stars we usually come across the limit T y m y a. Just to have an 
idea of the expression, we calculate the leading order contribution, at large 
chemical potentials due to the Higgs exchanging vertex, which comes out to 
be like, 



do* (it -> ii 
da (ti -> ii 



2aa „ l + v 
-in- h 



a 



E 3 



irv 2 E 1 — v Qirvm 3 



1 + in 



1-v 



(29) 



T y fj, could be important for the superdense collapsing stars like super- 
nova. It can be clearly seen that each term can have infrared divergences 
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in certain limit. We can play with the comparative limits to find physically 
interesting limits. The Z° vertex will not be any different. Only the coef- 
ficients of mass divergent terms like in would be different. This thing 
may give some help to address the stability of the superdense stars also. But, 
in such collapsing systems the densities are high enough to worry about the 
magnetic fields. We therefore, postpone this discussion for the moment and 
move over to the conclusion of our present results. 



4 Results and Discussions 

The results calculated in the previous section are mostly relevant from the 
point of view of cosmology. In the beginning of the universe, the temperature 
of the universe was close to the decoupling temperature, and the leptonic 
processes were taking place. Eq.27 shows the radiative corrections to the 
cross-section are simply proportional to T 2 which is the standard expected 
results in literature Ref. [18,19]. It is also worth mentioning that the above 
mentioned results are correct in the standard model. In the extended stan- 
dard model, we have to incorporate massive neutrino effects and several other 
contributions. 

The present experimental limit [20] on the standard model Higgs mass is 
m H y 107 GeV. In this situation, we cannot further simplify Eq.27, because 
we cannot ignore one term as compared to the other in bracket due to com- 
parable contributions. Just to have a very quick estimate of the magnitude 
of this contribution, we can assume the mass of Higgs to be 100 GeV, T 
=2MeV and the standard values of the gauge boson masses to use them in 
Eq.27. The radiative corrections from the background, then appears to be, 



da (ti -> it ) 



10~ iU (30) 



This contribution would be increased with increasing temperature but cannot 
go beyond a couple of orders higher, otherwise, the validity of the calcula- 
tional approximations would be questionable after T y a few tens of MeV. 
It happens because we deal with the exponential functions in the statistical 
distribution functions, and the approximations used for them are usually for 
a small range. 

Therefore, we reach the conclusion that FTD corrections to the scattering 
cross-sections are not so much significant in the early universe. But, the 
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energy density of the universe goes with the fourth power of these scattering 
cross-sections [3], so we still need to evaluate it to improve the precision in 
estimating the energy density of the early universe. On the other hand, the 
FTD corrections are not so insignificant after the universe was at a little 
lower temperature and the nucleosynthesis started taking place [22]. 
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